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Abstract 

We derive a non-linear supersymmetric cr-model for the transport of light 
(classical waves) through a disordered medium. We compare this model 
with the well-established non-linear cr-model for the transport of electrons 
(Schrodinger waves) and display similarities of and differences between both 
cases. We show that for weak disorder, both models are equivalent (have the 
same effective Lagrangean). This effective Lagrangean correctly reproduces 
the (different) Ward identities for Schrodinger waves and for classical waves. 

I. INTRODUCTION AND MOTIVATION 

In recent years, the transport of electrons and light through disordered media has been 
studied intensely, and many interesting effects have been observed and understood. Exam- 
ples are universal conductance fluctuations and weak localization for electrons, and speckle 
patterns and enhanced backscattering for light. A thorough discussion may be found in Refs. 

The universal tool to deal with these and other phenomena in the case of electrons 
has been Efetov's supersymmetric non-linear a-model (SUSIG) 0. This model success- 
fully describes not only the perturbative effects mentioned above but also non-perturbative 
features like localization. It correctly accounts for both transport properties and spectral 
fluctuations. Thus, it is fair to say that SUSIG embodies the essence of electronic proper- 
ties of disordered media. SUSIG has so far not been extended to the transmission of light 



through disordered media. In the present paper, we aim at filhng this gap. Our motivation 
for this work is the following. 

The transmision of light through disordered media is commomnly described in terms 
of the scalar wave equation rather than a variant of Maxwell's equations 0. The scalar 
wave equation differs in a fundamental way from the Schrodinger equation for electrons, 
see Sec. |l|. By the same token, the Ward identities for both equations differ substantially. 
In trying to extend SUSIG to the scalar wave equation, we probe the ability of the non- 
linear (T-model to provide a universal description of wave propagation in disordered systems 
described by different wave equations. It is of interest to see in which way the difference in 
wave equations is reflected in the effective Lagrangean of the resulting SUSIG. We will show 
that SUSIG does apply to the scalar wave equation, and that Efetov's effective Lagrangean 
is universal: It has the same form for Schrodinger waves and for scalar waves. Using the 
replica trick, John and Stephen |^ derived a non-linear cr-model for classical waves. This 
derivation was confined, however, to waves at fixed energy and thus bypassed the crucial 
issue of correlations between amplitudes at different energies. The latter play the central 
role in SUSIG. 

The paper is organized as follows. In section |I| we compare the wave equations for 
electrons and classical waves and derive the simplest variant of the Ward identities for 



both. Section III describes the derivation of non-linear a-model for classical waves in 



the supersymmetric formalism. We compare our result with the analogous expression in 



Efetov's work IBI. Section IV is devoted to the derivation of the Ward identities within a 



supersymmetric formalism for classical waves. Our conclusions are presented in section ^ 

II. WAVE EQUATIONS 

We put the mass m of the electron and Planck's constant h equal to unity. The 
Schrodinger equation for a noninteracting electron, 

(-lA + V^(r))0 = E0, (1) 



contains the random potential V{r) which describes impurity scattering. The propagation 
of hght is described by the classical wave equation 

(-A + k^6e{r))(f) = k^(j) . (2) 

Here, k = u/c is the wave number. We have decomposed the space-dependent dielectric 
constant e(r) = 1 — 6e{r) into a space-independent background term (which we put equal 
to 1) and a fluctuating part 5e(r). We assume that 6e is a Gausian random process with 
vanishing first moment and a second moment given by 

< 5e(ri)5e(r2) >= J^^i^i " ^2) , (3) 

where / is the elastic mean free path and d is the dimension of the system. We have written 
Eq. (Q) in complete analogy to the case of electrons. 

We compare Eqs. ([l|) and (H). Aside from a factor two, the quantity k"^ corresponds 
formally to the energy E. However, fc^ is always positive, in contradistinction to E. The 
main difference between Eqs. (|T]) and lies in the energy dependence of the random 
potential in the classical case. While V is independent of energy, the analogous term 6e is 
not. This difference is also reflected in different Ward identities which relate averaged one- 
and two-point functions. For electrons, the retarded (advanced) Green function 

Gf = iE^±z + A/2-V)-' (4) 

is taken at energy E^ + z = E + z + iri {E~ — z = E — z — if]), respectively. We immediately 
find 

<Gt>-<G: >= -2{z + ir,) < GjG; > , (5) 

where the angular brackets stand for the ensemble average. In complete analogy, we define 
the Green functions for classical waves by 

(6) 
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The Ward identity reads 

< Gt > - < G- >= -2 + irij < GtG- > +Ak^ < G+^eG; > . (7) 

Because of the frequency dependence of the impurity term in the classical case, these two 
Ward identities differ in form. They actually also indicate different conservation laws: Parti- 
cle conservation for electrons, and energy conservation for classical waves. The Ward identity 
for the classical case will serve as a check of our supersymmetric formalism: In section |^ 
we derive it from the non-linear a-model. 



III. NON LINEAR SIGMA MODEL 

We derive the non-linear cr-model for the simplest non-trivial case, the ensemble average 
of a product of an advanced and a retarded Green function. We use the notations and 
definitions of ref. f^. The advanced and retarded Green functions can be written as integrals 
over supervectors 

G±(yi,y2,fc') = T^/i?[^]^.(yi)^I(y2)exp[/:(vi>)] , (8) 

where we have omitted the index c for the Green's function of classical waves, and the 
Lagrangean is given by 

C = ltJ d^y [vl>t(y) (± + A - 5t{y)e) + zr/) v]>(y)] . (9) 

The quantities ^'(x) are supervectors defined by 

v[/(x)t = (^i(x), ^2(x), -x(x) , x*(x)) . (10) 

The quantities S are ordinary real integration variables, and the x's anticommute. We in- 
troduce a source term J(y) =diag(j(y), 0, 0, 0) in graded space and introduce the generating 
functional 

Z^{k\ J) = j exp U + d'y{^Ky)J{y)^{y))] . (11) 



This functional generates the Green function at point yi = y2, which is sufficient, because in 
the present section we are interested only in the effective action. In section IV we show how 
to generate the Green function with different space point arguments. The Green function is 
given as functional derivative of the generating functional with respect to J at J = 0, 



dj{y) 



(12) 



We use this expression to calculate the average of the product of a retarded and an advanced 
Green function taken at different frequencies, < G^G^ > (the two-point function). This 
quantity plays an important role in describing average properties of random systems, such 
as the level-level correlation function, the distribution function of the transmission, etc. It 
serves as an example. Except for the dimension of the Q-matrices appearing below, and 
except for the dependence on additional frequency variables, the average 2/c-point function 
for any positive integer k is governed by an effective Lagrangean of the same type. 
The generating functional Z for the two-point function is given by 



Z{P,Ak^,J) 



P[^]exp 



vl/t(y)LV2 (kl + + A - 5e(y) (k^ + 



+ J(y))Li/2^(y) 



(13) 



where L =diag(l, 1, 1, 1, —1, —1, —1, —1), \1/ are supervectors with 8 components, and J is 
an 8 X 8 matrix. All quantities are given in "advanced-retarded" notation (see Ref. ||^). 
Averaging over the Gaussian distribution of 6e, we obtain the Lagrangean 

£ = 1^ / d'y vl/t(y)LV2 (^kl + ^L + A + ^vL^ ^(y)^^/^- 



2k 



V(y)LV2(l + ||L)LV^vl;(y)y 



(14) 



Using the Hubbard-Stratonovich transformation in the usual way and integrating over the 
vectors we obtain the following form of the generating functional. 



VQ exp 



Txu 1 / AA;^ 

-— trgg2 + -trglog [kl + + A + zr^L + J(y)- 
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1 / Ak 



(15) 



Here v is the density of states per unit of and per unit of volume, and r = fcQ^^£/(27r^z/) 
formally corresponds to Efetov's mean free time 0. We have introduced these quantities in 
Eq. (|T5|) in order to facilitate the direct comparison to Efetov's expression for electrons. The 
term (l/2r)QL(A/c^/2A;Q) is due to the frequency dependence of the "scattering potential" 
k'^5e in Eq. (|^). Comparing Eq. ( |I5| ) with the corresponding expression in Efetov's work 
1^, we identify (modulo factors of two) e^k"^ with the sum energy and e^Ak"^ with the energy 
difference and find that the two expressions differ by the term (l/2r)QL(A/c^/2/cQ). 

To evaluate Eq. (|l^), we use the saddle-point approximation. This is justified if r <^ p, 
the mean level density. Varying the Lagrangean in Eq. ([T5|) with respect to Q and neglecting 
terms proportional to Ak"^ and source terms, we obtain the standard saddle-point equation 

(16) 



Q = — tr 
vrp 



kl + A-^Q 



This is the same equation as in the case of electrons. As in that case, the condition koi 3> 1 
(weak disorder) yields Q = iL as a solution of the saddle-point equation. The weak disorder 
condition also implies, however, that the term (l/2r)(5L(AA;^/2/cQ) in Eq. (|15D can be 
neglected. This is the case for sufficiently large ko. Then, there is no difference between the 
non-linear a -models for Schrddinger waves and for classical waves. This statement is the 
central result of our work. It obviously extends to the generating functionals of all higher 
correlation functions and, thus, applies universally. 

The actual differences between the two theories are due to the different forms of the 
source terms. In the next Section, we show this in the case of the Ward identities. 



IV. WARD IDENTITY 

In Appendix G of Ref. [Q, it was shown how a Ward identity can be derived in the 
context of SUSIG. We use that method to check the Ward identity, Eq. (|^, for classical 
waves, using essentially the generating functional derived in the previous Section. With 
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slight modifications, our calculation also applies to the case of electrons. We first show how 
the new source terms emerge, when we introduce a new generating functional for the r.h.s. 
of Eq. d^). We use the coordinate representation 



< r|G+(l -(5e)G-|r' >= J rf^xG+(x, r)(l - 5e(x))G-(r', x) . 



(17) 



The generating functional 



Zi(A;^ AA;2, Ji) = 

y /dW^I(y)^i(y,y')^i(y' 



^i(y') + 



(18) 



produces the retarded Green's function on the r.h.s. of Eq. ([TTI): 

G^(^'r) = -^, (19) 

where the source term Ji = diag(ji, 0, 0, 0) is put equal to zero after taking the deriva- 
tive, while is a supervector with 4 components. We also introduce another generating 
functional 



P[vl>2]exp[-^|rf^yvl>t(y) 



f/VV^ky)(l - 5e(y))J(y,y')^2(y')] • 



(20) 



Then, immediately 



;i-5e(x))G-(r',x) 



dZo 



(9j2(x, r') 



(21) 



Taking the product of Eqs. (p!8|,|20|) , we obtain a generating functional 



V[^] exp 



(22) 



where the action is given by 



Jiy,y')-Seiy')Jiyy)p/'^iy') , 



(23) 



with J = diag(0, 0, 0, 0, j2, 0, 0, 0), J = diag(ji, 0, 0, 0, j2, 0, 0, 0). The second partial deriva- 
tives of Zf produce the integrands on the r.h.s. of Eq. (|^). The additional source term 
J(y, y') represents the important difference to the electron case. Averaging of the term 
containing the random part of the dielectric constant 6e{y) leads to 

(24) 



exp 



where 



^(y,y') 



%-y') + ^(y,y') 



(25) 



To perform the Hubbard-Stratonovich transformation we introduce a supervector 



i^(y) = /rfV^(y,y')L^/'^(yO 



(26) 



after which we can rewrite the expression in Eq. (^4|) as 

1 



exp 



ciV[^^(y)L^/'F(y)]^ 



(27) 



Using the Hubbard-Stratonovich transformation and keeping only diffusive modes, we obtain 
the average of Eq. (p^) 



Zf = Vl"^] exp 



= V[Q]exp[-C,iQ)] , (28) 



dV?^trgg2 _ i logDetg5(g) 



where Detg means determinant over real and graded spaces and we define a matrix 



B{Q) 



H L + A + ir]L 



^(y-y') + ^(y,y') 



2Tk^ 



Q(y)A(y,y') 



(29) 



For maximum compactness we allow J and J to be general symmetric 8x8 matrices. 
This is permissible, because we never use the particular form of the source term in the 
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derivation in the previous section. The saddle-point equation for this action is the same 
as before. Following the formalism developed in Ref. |^ we apply the transformation Q 
(1 + 5T)^^Q{1 + 5T) (we preserve the notation of ref. |^), changing the action in Eq. ( |28D 
into 



c, + ]p:igB-\Q) 



ap 

^T. ^U{y - y') 



5T,(^L + zr^L)5(y-y') + J(y,y') 



^~(y, y') 



(30) 



where Trg is the trace in both, real and graded, spaces. A transformation of integration 
variables leaves Zj invariant. Therefore, terms linear in 6T in the expression of Zf must 
vanish, which leads to the following equation 

6T,J 



V[Q] exp{-C,)TrgB-\Q) 
+ J V[Q]exp{-C,)TTgB-\Q) 



AP 



2Tkl 



V[Q] exp(-A)TrgS"i(g)g 



AP 

5T,((^L5(y-y') + J(y,y')) 



(31) 



We will consider each of these terms in detail. The first term is 



P[Q]exp(-A)Trg[i?-i(y,y')]"^ 



5T,J 



y.y 



J2j J dVV/ v[Q] exp(-A)5-i(g) 

af3 



5T,J 



Pa 



kk' 



,dZj{klAk^J) 



5T(1,2)J(2,1) - J{1,2)ST{2,1] 



kk' 



y ,y 



kk' 



+ E/ d'y'd'y 



dZf{klAe,J) 



d4Tiy',y) 



5T(2,1)J(1,2) - J(2, l)5r(l,2) 



kk' 



+ ... 



(32) 



y .y 



where we use block notation as in Ref. 0. The dots represent terms containing J(l, 1) 
and J(2, 2) which do not contribute to the final result. Using the explicit expressions 

'^^fefc' = "^^kko^k'k'^'' ^^^kk' = -'^^i.i.'h'ka, the first term can be written as 



'kk'^"k'ko, 



V^^'^'^ <\V(y',y) V^^'^) 
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34l:"(y',y) 



kka 

.2 A U2 



'i://«1f^4:?(y',y, 

_.x://.vvfi^.s:'(y',y).... 



I'-i OA 

Now taking the derivative with respect to Jf^' (x, x'), we finallly obtain 



2i 



dZf{klAk^,J) 



J,J=0 I'^o 



5i , ' - 2i 



dZf{klAe,J) 



aj(?f(x,x') 



J,J=0 



We consider now the second term in Eq. (^Tj) 



exp(-£i)1VgB-'(g)|ir, (— L + «,L)*(y - y')| = 

+2(^ + E / / AciV^?4S^^^5T.„(2, l)5(y - y') 



kk' 



,Ae , f dZf(klAk\J) .AP ^ /• , dZfikl,Ae,J) 



The derivative with respect to Jii'^'*(x, x') leads to 



^Ak^ , / f d^Zf(klAk^,J) 

-2i( + ir])[ cf^y — °' I ' 

2 '^9j(;;J)(y,y)9Jli'^)(x,xO 



/A 



d^Zf{klAk\J) 



9j|;^;^)(y,y)9J}r(x,xO 



J,J=0 



J,J=0 



+ 



We expand the remaining term in Eq. (|31|) 
1 



2Tkl 



,AP 



V\Q] exp(-A)Trg5-i(g)g[5r, (^L5(y - y') + J(y, y')] 



.,dZf{klAk\J) 



6T. 



Ak 



2 "I 

-L 



zAfc^ f d\??Mi^^ _ f^^ 9Z-j{klAk\J) 
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The term proportional to J is omitted, because it does not contribute to the final result. 
Taking the derivative with respect to \ we find 



d^Yjikl Ak\ J) 



J,J=0 



J,J=0 



(38) 



In conclusion, the requirement that the term linear in ST in the expansion of Zf vanishes, 
entails the following equation : 



dZfiklAk^J) 



dZf{klAk\j) 



2 +^^V^^,j(?^)(y,y),ja^2)r_^ 
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J,J=0 



X,X' 



J,J=0 



-2i{ 



r ^ ^^ZA^^Ak^Jl_ 



J,J=0 



d''Zf{kiAe,j) 



+Ak^ I d^y-^ 



dJkok',^'''\y,y)ji'{'\ 

d^Tf{ki,Ae,j) 



X, X 



;y,y)jS''Hx,x' 



J,J=0 



J,J=0 



(39) 



Replacing Zj everywhere by Z, using definition of Z and putting k^ = k'Q = Iwe immediately 
obtain a Ward identity in the form of Eq. (|^). Indeed, Eq. (Q) corresponds to the l.h.s. of 



Eq. (0), whereas Eqs. 



lead to the r.h.s. of Eq. (|^). 



V. CONCLUSIONS 

We have derived a non-linear supersymmetric a-model for classical scalar waves. We 
have shown that in the weak disorder limit {k^l ^ 1), the effective Lgrangean of this model 
is identical to the one for electrons. In this limit, the main difference between the wave equa- 
tions for classical and Schrodinger waves, the frequency dependence of the random potential. 
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does not lead to different wave behavior. We have also shown that the Ward identities for 
classical and for Schrodinger waves are both fulfilled by the same effective Lagrangean. This 
is due to the different source terms. Outside the regime of weak disorder, i.e. at low fre- 
quencies, the frequeny dependence of the disorder potential for classical waves suppresses 
disorder effects altogether. This does not happen for electrons. This low-frequency domain 
is not accessible to the non-linear cr-model. 
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